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We study QCD with two staggered Dirac fermions both in the fundamental (QCD) and the adjoint 
representation (aQCD) near the chiral transition. The aim is to find the universality class of 
the chiral transition and to verify Goldstone effects below the transition. We investigate aQCD, 
because in that theory the deconfinement and the chiral transitions occur at different temperatures 
T c i < T c . Here, we show that the scaling behaviour of the chiral condensate in the vicinity of 
j3 c is in full agreeement with that of the 3d 0(2) universality class. In the region T c / < T < T c 
we confirm the quark mass dependence of the chiral condensate which is expected due to the 
existence of Goldstone modes like in 3d 0(N) spin models. For fundamental QCD we use the 
p4-action. Here, we find Goldstone effects below T c like in aQCD and the 3d 0(N) spin models, 
however no 0(2) /0(A) scaling near the chiral transition point. The result for QCD may be a 
consequence of the coincidence of the deconfinement transition with the chiral transition. 
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1. Introduction 

We study QCD with two staggered Dirac fermions both in the fundamental (QCD) and the 
adjoint representation (aQCD) near the chiral transition. The intention is to find the universality 
class of the transition and to verify Goldstone effects below the transition. For ordinary QCD with 
staggered fermions the prediction of the three-dimensional 0(2) universality class (the 0(4) class 
in the continuum theory) could up to now not be confirmed. We investigate in addition aQCD, 
because in that theory the deconfinement and the chiral transitions occur at different temperatures 
with Td < T c (that is /3d < j3 c ) [fi]|. The chiral transition can therefore be studied without interfer- 
ence. The comparison of QCD data with the critical behaviour of 0(N) spin models requires the 
following identifications of QCD variables to 0(N) variables: the chiral condensate corresponds to 
the magnetization M and the quark mass (m q a) to the magnetic field H. Instead of the temperature 
we use /3 = 6/g 2 . 

2. QCD with adjoint fermions (aQCD) 

The action of aQCD which we use is ^ 

S = f3S G (U^) + Y,V* M ( u{ * ) ^yVy • ^ 

Here, the gluon part Sg(U^) is the usual Wilson one-plaquette action, but the fermions are in the 
8-dimensional adjoint representation of color SU(3). The standard staggered fermion matrix M 
depends correspondingly on instead of U^ 3 \ The links are real because 
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(2.2) 



The fermion action does not break Z(3) center symmetry and the Polyakov loop L3 is therefore 
order parameter for the deconfinement transition. In the continuum aQCD contains an SU (2Nf) 
chiral symmetry which breaks to SO(2Nf) for T < T c . This is because here Nf Dirac fermions 
correspond to 2Nf Majorana fermions. For Nf = 2 we have SU (4)-symmetry which breaks to 
SO (4). The corresponding continuum transition has been studied in Ref. [§] with renormalization- 
group methods. On the lattice an (?(2)-symmetry remains for staggered fermions. 

Our simulations [^] were done on N a x 4 lattices with N a = 8, 12 and 16 and a fixed length 
T = 0.25 of the trajectories. We used 900-2000 trajectories for measuring the chiral condensate 
(YY), the Polyakov loop L3 and the disconnected part %dis of the susceptibility. 

In Ref. ^ the deconfinement transition point was found at = 5.236(3) - the same value as 
in QCD(\), but here it is of first order. The usual strategy to locate the chiral transition point j8 c is 
to extrapolate the line of peak positions j8 pc of the susceptibility to m q a = with 

P pc = Pc + c(m q a) l / A , (2.3) 

where A = jS m 5 is a product of critical exponents. On the 8 3 x 4 lattice we find for m q a = 
0.005,0.01,0.02 the values ji pc = 5.73(8), 5.74(6), 5.75(10). With A from 0(2) we obtain j8 c = 
5.7(2). A better method [^] is to expand the scaling ansatz {yy) = ml ^ S f(z) a t at small |z|. Here, 
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Figure 1: The chiral condensate (iff iff) of aQCD as a function of (m q a) 1 / 2 for all j3-values between 5.3 
(highest values) and 5.9 (lowest values) from jV^x4 lattices with N a = 8 (circles), 12 (triangles) and 16 
(diamonds). The black lines are fits with ansatz (2.5), the filled circles denote the extrapolations to m q a = 0. 

m = m q a/(m q a)o, j8 r = (j8 — j8 c ) / j8o are the reduced field and temperature and z = j3 r m~ l / A . The 
result is 

(xff ¥ ) = ('n q a) 1/S {d c + d^-p c )(m q a)- l / A + ...} . (2.4) 

Fits to only the first term at fixed /3 are best for /3 € [5.6,5.7]. If 0(2) exponents and two terms are 
used one finds a unique zero of the parameter d^(P — /3 C ) as a function of /3 at the rather precise 
value p c = 5.624(2). 

Due to the existence of massless Goldstone modes for all T < T c the susceptibility %l of 3d 
0(N) spin models diverges for H — > as Xl ~ H~ x l 2 . If aQCD behaves effectively as such a 
model, a similar divergence is expected. In terms of the chiral condensate we must have then 

(yy) (0 , m q a) = {xj/yf) (j3 , 0) + a (J3) {m q a) l l 2 + c 2 (P)(m q a) + ... (2.5) 

Corresponding fits are shown in Fig. [I] for j8 = 5.3 to /3 = 5.9 with A/3 = 0.05. The blue line 
separates fits above and below /3 C . We see that the model expectations are met well by the data. 

We have performed explicit scaling tests for the data from the 8 3 x 4 lattice (they coincide 
with the results from the N a = 12 and 16 lattices). In the left part of Fig. ^| the data for (l/Zi/^ra -1 / 5 
are shown as a function of the scaling variable z = fi r /m l l A using 0(2) exponents. The normaliza- 
tions (m q a)o and j8o have been determined from the behaviour at /3 C ( |2.4[ ), and the extrapolations 
(\j/Y)(P,0) obtained from the Goldstone effect [ph. In addition we show the 0(2) scaling function. 
We observe that the data scale very well with 0(2) exponents in the small \z\ -region and agree there 
with the 0(2) scaling function. For decreasing z < outside the shown range the data for different 
/3 start to deviate and exhibit corrections-to-scaling as in the original 0(2) model. In the right part 
of Fig. § we show the best scaling result for the unnormalized variables obtained by varying the 
exponents around 8 « 4.0, v « 1.1, the expected continuum class values Here, /3 f = 5.64. We 
see from Fig. || that the data still spread considerably, a scaling function is not known up to now. 
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Figure 2: Scaling test in aQCD. Left part: {xj/^/m 1 ^ 5 for 0(2) parameters in the vicinity of the critical 
point. The black line is the 0(2) scaling function. Right part: the unnormalized variables with 8 = 4.4, v = 
1.0 from Ref. [§], The data for fixed j3 -values are connected by straight lines to guide the eye. 



3. QCD with fundamental fermions 

Since the predictions for the universality class of ordinary QCD could not be confirmed with 
the standard staggered action we use here the /?4-action[f|], which improves the cut-off dependence, 
rotational invariance and flavour symmetry. The thermodynamics of two flavour QCD has been 
investigated with this action in Ref. In particular, the chiral transition point was estimated to 
j8 c = 3.48(3) by extrapolation of the pseudocritical points fi pc (m q a) to m q a = 0. 

We have extended the work of Karsch et al.[]|] on N„ x 4 lattices with N a = 8, 12, 16 at 
the couplings j8 = 2.8,3.0,3.2, ,3.4,3.48 and 3.50, that is for j8 < j8 c . All our simulations were 
done with the MILC code using the R-algorithm with a mass-dependent stepsize 8z(m q a) = 
min{0.4m g a,0. 1} and trajectory length z = 1. We produced 1000-2000 trajectories for masses 
m q a 6 [0.025,0.5] and somewhat less for m q a = 0.01. 

Like in aQCD we have tested the mass dependence of the chiral condensate at fixed j8 < j8 c . 
In Fig. H we show the data from j8 = 2.8 to j8 = 3.48 as a function of (m q a) l l 2 . We have fitted the 
data in the range m^a<0.15 to the ansatz 

(W>(/3,m,a) = <W>(JS,0) +c 1 ( j 8)K«) 1 / 2 . (3.1) 

Again, the behaviour expected due to the Goldstone effect is confirmed by the data. 

In Fig. ^ we investigate the scaling form = d c (m q a) l ^ s for the critical point in the close 
neighbourhood of j8 c . We have plotted the data with 8 (2) = 4.78 in the left part of Fig. |] at j8 = 3.48 
(solid lines) and j8 = 3.5 (dotted lines). Obviously, the curves are no straight lines going through 
the origin. Also, the corresponding attempt to obtain a scaling function with 0(2) exponents fails 
as can be seen in the left part of Fig. |[ Here, /3 C = 3.51 was used, the scales are not normalized. 
Other values for j8 c shift the results slightly, but do not lead to scaling. The use of 0(4) exponents 
does not improve the result. 
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Figure 3: The chiral condensate (iff iff) as a function of (m q a) 1 / 2 for all f3 -values between 2.8 (highest values) 
and 3.48 (lowest values) from N% X 4 lattices with N a = 8 (circles), 12 (triangles) and 16 (diamonds). The 
black lines show the fit results and the filled circles the estimates for the chiral condensate at m a a = 0. 



In order to achieve scaling of the data we have fitted the ansatz (iffiff) = d c (m q d) x l^ at /3 = 3.50 
with a free amplitude d c and a free exponent 8 for m q a € [0.01,0.1] and obtain 8 = 2.3. The fit is 
shown in the right part of Fig. |[ Likewise, we estimated the magnetic exponent j8 m to 0.6 from the 
values of (iffiff) at m q a = which we gained from the Goldstone extrapolations. In the right part of 
Fig. |5| we show the data for these exponents. Here, (iffiff) / (m q a) 1 ^ 5 is normalized to 1 at j8 c . Still, 
scaling is not perfect, but certainly much better than for 0(N) exponents. 
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Figure 4: Test of the scaling ansatz for (iffiff) at the critical point for 0(2) (left part) at j3 = 3.48 (solid lines) 
and 3.5 (dotted lines). In the right part we show a free fit to the scaling ansatz at j3 =3.5 (dashed line). Here, 
the exponent is 8 = 2.3. 
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4. Summary 

In aQCD the deconfinement transition is first order and occurs below the second order chiral 
transition. The latter is located at /3 C = 5.624(2). The scaling behaviour of (yy) in the vicinity 
of j6 f is in full agreement with the 3d 0(2) universality class. The lattice data do not yet show 
scaling with exponents from the proposed continuum class [Q]. In the region between the two phase 
transitions the quark mass dependence of the chiral condensate is as expected due to the existence 
of Goldstone modes like in 3d 0(N) spin models. 

For QCD we find that the chiral condensate exhibits Goldstone effects below the chiral tran- 
sition point as in 3d 0(N) spin models and like in aQCD. The transition is of second order or a 
crossover, there is no sign of a first order behaviour. In the vicinity of j3 c the chiral condensate does 
not show the scaling behaviour of the 3d 0(2) or 0(4) universality classes, at best it is that of a 
class with exponents 8 = 2.3 and j8 m = 0.6. This result could be a consequence of the coincidence 
of the deconfinement and the chiral transitions. 
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